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We discuss the topology of the wavefunctions of two-dimensional time-reversal symmetric 
■ superconductors. We consider (a) the planar state, (b) a system with broken up-down reflection 

JzL ■ symmetry, and (c) a system with general spin-orbit interaction. We show explicitly how the relative 

^ | sign of the order parameter on the two Fermi surfaces affect this topology, and clarify the meaning 

^v^j ■ of the Z2 classification for these topological states. 
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I. INTRODUCTION 

' Topological states of matter have received much recent attention^ Best known examples are the two-dimensional 
integer quantum Hall states.— These states can be characterized by the Chern numbers associated with the electronic 
wavefunctionsi 3 - At the boundary of the sample, the non-trivial topological structure of these wavefunctions necessarily 
implies the existence of chiral edge states. These states can carry charge (Hall) currents even though the bulk is gapped. 
. Chern numbers can be finite in systems only with broken time reversal symmetry, which in this case is provided by 
the external magnetic field. Recently, there is much interest in time-reversal symmetric topological insulators4£& 7 -'£ 

S These insulators cannot be classified by the Chern numbers, yet they are topologically distinct from ordinary insulators 
we learnt in standard solid state textbooks. The role of the Z2 topological number in these systems has been elucidated 
1-^ , recently, especially by Kane and his collaborators^ 5 -! 6 -. There are also surface bound states at the edges of a topological 
insulator. For a two-dimensional topological insulator, there exist a pair of one-dimensional edge states, related by 
time-reversal symmetry, at the boundary of the system. In the absence of spin-orbit interaction, these edge states 
can be regarded as two one-dimensionally dispersing states of opposite spins traveling with opposite velocities. These 
pair of states can carry spin current, even though the bulk states are gapped. These topologically insulating states 
£SJ . are often also referred to as spin-Hall insulators^ 

' There are directly analogies of topological states in superconductors. Here, while the system possesses an extra 
mechanism of transport via motion of Cooper pairs, the quasiparticle wavefunctions can have topological properties 
On . imposed upon them by the order parameter similar to those discussed above. Neither is one restricted to a periodic 
system; the discussion can be done also in the continuum case. A direct analogy to the integer quantum Hall states are 
spinless weak-coupling 11 superconductors where the order parameter are (k x + ik y ) n , where k x>y are the components 



of the wavevector k and n's are (odd) integer o 12 ' 13 ' 14 ' 15 . For the continuum, there is no limit to the magnitude of k. 



The Chern number is n. For n = 1, we have the familiar two-dimensional p wave superconductor which possesses 
a chiral gapless state at the edge of the sampl o 13 ' 15 . The non-trivial topology of this state is also reflected in the 
bound state spectra in a singly quantized vorte x 14 ' 16 . In this case, one finds a Majorana mode at exactly zero energy. 
Majorana vortex bound states and edge modes are intimately related to each othe r 15 ' 17 . The direct analogy of a 
time-reversal symmetric topological insulator in two-dimension is the state often referred to as the "planar state" in 
the 3 He literature 1 ^. Physically, this state is a " equal spin pairing state" where the up and down spins pair only 
among themselves (but not with each other), with the down (up) spins paired into a (k x + ik y ) ( (k x — ik y )) state. 
This state is obviously time-reversal symmetric. Since the up and down spins are decoupled, it is also clear that the 
boundary of such a two-dimensional system would have one edge state for each spin species, propagating in opposite 
direction s 19 ' 20 (see also 24 -). A singly quantized vortex in the planar state has two Majorana modes, one for each spin 
species. In contrast, an s-wave superconductor is topologically trivial. It has no topologically required edge states at 
the sample boundaries. No Majorana mode exists in its vortex core, and the lowest energy bound state has a finite 
energy^. 

What is the fate of this superconducting state when the spins are no longer decoupled? Besides basic theoretical 
curiosity, this question has now an added significance due to the study of non-centrosymmetric superconductors. 
These systems include three dimensional systems such as CePt^Si and related compound s 23 ' 24 , or two-dimensional 
systems at the interface between two different materials 25 . Consider here in particular a two-dimensional system 
without up-down reflection symmetry, though still with rotation symmetry about the normal to the plane. In this 
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system, the planar state is indistinguishable from the s-wave pairing state by symmetry, so the two order parameters 
can coexist^. However, an s-wave state is topologically trivial, so there must be a non-trivial interplay between the 
two order parameters when they are simultaneously present. In an earlier paper—, we investigated the presence or 
absence of zero energy modes at the vortex core for such a state. There we find that the zero energy modes exist 
only if the p-wave order parameter is larger than the s-wave order parameter at the Fermi surface. In addition, we 
have included a Rashba energy term in the single particle Hamiltonian, which spin-splits the Fermi surface into two 
according to helicity. We concluded in this case that the existence or absence of the zero energy states is determined 
by the relative sign of the pairing order parameter on these two Fermi surfaces^ 

The edge states at the boundary of a two dimensional non-centrosymmetric superconductor have also been studied 
ir*2&. There the authors have found the conditions on existence or absence of edge states consistent with what we 
have stated above, (see also^ for further discussions) 

In this paper, motivated by the above, we try to understand the results of re f 2 ^ 28 by directly examining the topology 
of the quasiparticle wavefunctions in the bulk. A fundamental concept in topological states is "obstruction", that is, 
the inability to construct a wavefunction smooth in some parameter space (here momentum space) satisfying certain 
basic criteria. For the integer quantum Hall state, this has been illustrated explicitly by Kohmoto^S. There he showed 
that, for electrons in a lattice and in the presence of a magnetic field (with rational flux quanta per unit cell), one 
cannot have a Bloch wavefunction that is smoothly defined throughout the entire (magnetic) Brillouin zone. Dividing 
the momentum space into two regions, he then showed that the transformation function relating these two regions 
(also known as the "transition function") reflects the finite Chern number of the state. Similar consideration has 
been utilized by Fu and Kane (Appendix A of 5 ): an insulator is topologically non-trivial if one cannot construct a 
pseudospin basis smooth within the Brillouin zone (more precise statements below). A transformation function is 
then used to define the Z2 topological number. We shall follow the same route here for our superconducting states, 
with k however which can extend to infinity. We shall see directly how the relative sign of the order parameter on 
the Fermi surfaces determines whether it is possible or not for the wavefunction to be smooth in k space. This thus 
in turn determines the Z2 topological classification of the superconducting state. 

While preparing the present manuscript, a preprint^! on the topology of time-reversal symmetric weak-pairing 
superconductors in general dimensions appears, employing a method that is very different from ours. These authors 
started from the quasiparticle Hamiltonian rather than the wavefunctions. For three-dimensions, they developed their 
classification by employing a topological invariant introduced earlier 10 . Their classification in two-dimensions is then 
obtained by performing a dimensional reduction from three-dimensions. For multi-Fermi surfaces, they showed that 
the product of the signs of the order parameter on Fermi surfaces enclosing time-reversal invariant points determines 
which Z2 class the superconductor would belong. For a single pair of Fermi surfaces enclosing k = 0, their result 
reduces to what we have stated above. The arguments in Rcf 31 , though general, are in a rather abstract form. For 
two-dimensions, they also rely on an artificial extension of the physical system to an extra dimension. We hope that 
some readers would find our explicit consideration here of obstruction and transformation function within the physical 
system more easily understandable and informative. 

Our paper is organized as follows. For clearer explanations, we first start with the pure planar state, then introduce 
s-wave mixing and Rashba interaction. We shall explicitly construct the Z2 invariant. Finally, we generalize to a 
general two-dimensional time-reversal symmetric state. The Appendix contains some technical details. 

II. OBSTRUCTION AND THE TRANSFORMATION MATRIX 

A. Planar phase 

We begin by considering a superconductor (superfluid) in the planar state. The Hamiltonian is given by H = 
Hk + H pa i r where the kinetic energy 

k,oc 

and the pairing term H pair — > H*^* , with 

k 
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fc, 4>£ are the magnitudes and azimuthal angles of fc, a =1, | the spins, a is the corresponding annihilation operator, 
and h.c. denotes the Hermitian conjugate, m is the mass and fj, is the chemical potential. A p (fc), which depends 
only on the magnitude of k, is chosen real. In eq ([2]), the up (down) spins have angular momentum =pl. To have a 
Hamiltonian smooth in k space, we shall require that A p (k) — ► as fc — > 0. We shall adopt the following convention 

for the time reversal operator 0: Oat _1 = a - , and correspondingly Sat _1 = — - . The Hamiltonian 

fc,T -k,l * & J k,l -fc,T 

is time reversal invariant, for example 0ie _z< ^at - _1 = ie^Sat a -. . In principle we can also choose an 

ft, J — fc,f kA —k,l 

Hp l ^ m which is time-reversal invariant under only after an additional gauge transformation, but this would not 

add any new physics and we avoid such complications. We could also have chosen a different Hp^** 1 such as 

Ap(k)e~ l ^kat t — A p (fc)e i< ^Eat a' =• , but this would amount to replacing </>£ in eq ((!]) by <fe — § (that is, a 

relative spin-orbit rotation. This state corresponds to the one used ir>2£). Thus eq ([2]) is the most general planar state. 
The i factors in this equation are chosen for convenient comparison with the case with finite spin-orbit interactions 
(subsection B) below. 

This Hamiltonian can be easily diagonalized by Bogoliubov transformation, since the two spin species are decoupled. 



The energies are given by ±E(k) where E(k) = y£, 2 (k) + A 2 (fc), each doubly degenerate due to spin. Here £(fc) = 

J^- — [i. The unoccupied positive energy states 3 ^ |^^(fc) >, \^>^{k) > have wavefunctions^ 3 . given by, in the standard 
notation in spin and particle-hole space, 



*f (k) = N 



A 



1 



-iA p (k)e^/(E(k)+ak)) 




(3) 



9f (fc) = N J 



\ — iA p (fc)e 




1 


i/(E(k) 



(4) 



Here N A is a normalization constant, which we choose to be real and thus N A = ((E(k) + £(k))/2E(k)) 1 ' 2 . For a 
given fc, \^ A (k) > and \^! A (k) > form a complete set for the positive energy states of the Hamiltonian H. In writing 
eqs ([3]) and (|4]) we have required that they satisfy the criteria for a good pseudospin pair as stated in Appendix of 5 
(their eq (Al)), that is , 



and correspondingly 



0|*i(fc) >= |* a (-A) > 



e|* 2 (fc) >=-|* 1 (-fc) > 



(5) 



(6) 



when 1 ->| and 2 

We shall limit ourselves to the case where A p (fc) grows slower than fc 2 as k 
in this limit. Then > and \^ A (k) > are both well-defined as fc - 

are ill-defined at the origin fc — >■ if the chemical potential \x is positive 
(E(k) - £(fc))/A p (fc) -> oo, since £(fc) -> -ft < and A p (fc) 



0. Hence in this limit 36 
(fc) -» (0,0,-ie^,0) 



> oo, so that A p (k)/(E(k) + Z(k)) -> 
oo2^. However, these wave functions 
As fc -»■ 0, A p (k)/(E{k) + £(fc)) = 



*f(fc)->(0,0,0,-ie-*s) 



obviously ambiguous at fc = 0. 

We can however make the alternative choice, 
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*f (Je) = N 1 



/iA p (k)e-^/(E(k)-ak)) 

1 

V o 





iA p (k)e^/(E(k) - £(k)) 



(7) 



satisfying the requirements eqs © and ©, Here N B = ((E(k) — £(k))/2E(k)) 1 ^ 2 . These wavefunctions are well- 
defined at the origin. As k -> 0, A p (k)/(E(k) - £(fc)) -> 0, so 

*f(*)-> (0,0,1,0) 

*f(£)-> (0,0,0,1) 

However, they are ill-defined as fc — > oo. In this limit, A p (k)/(E(k) — £(&)) = (E(k) + £(fc)) / A p (k) — > oo so 

*f (fc) -» (ie-^s, 0,0,0) 

*f (A) -» (0,ie < *«,0,0) 

and have values which are 4>t dependent.— 

Thus (k) > forms a good pseudospin basis for any k excluding the origin, while |*^(fc) > is good for any k 

except at infinity. To have a good pseudospin basis for all k, one must divide the k space into two parts, say k > K 
and < k < K for some K, and use respectively \^^(k) > and |$^(As) > in these two regions. For our specific 
example, K can assume any finite non-zero value. At the boundary between these two regions (here k — K), we 
consider the unitary transformation matrix t AB between these two basis (see Appendix A o^, though the details of 
our analysis is different from this reference; see als o 37 ' 38 ), with matrix elements defined by 

^(^)=<^(^)|*f(0)> fc=JC (9) 

which is now a function of the angle <f> (with K a parameter). The requirements eq ([5]) ((6]) impose a non-trivial 
constraint on this matrix. We have 

(l*i(0 + tt) >, \M4> + »r) >) k = (|e*i(^) >, |e* 2 ($ >) fc ( \ J ) (io) 



for both A and B and 



We therefore find 



t AB ( 



, x f <*^ + tt)|\ (|^f + >,i*f C^ + tt) >) fe 



t AB (0+-) = (? v) tAB ^(-°l I) w 

Since t AB is unitary, we can write t AB = e tx ^'G((j>) where G is an element in SU(2) (detG = 1)^ and G(</>) is 
required to be smooth. Hence we have 

e 4xW+7r) G(</> + tt) = e _i * w) G(0) (13) 

and so there are two possibilities, either 

e *x(4>+*) = e -ixW + tt) = G(0) (14) 
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e ix(<t>+*) = _ e -ixW G((f> + tt) = -G((f>) (15) 

That is, if we consider the map from (j> £ [0,7r] to G((j)), the SU(2) part of t AB , it must be either periodic (eq 
(fT"4]) ) or antiperiodic (eq (jT5j) ) . We note that, if this map is antiperiodic, there is no smooth deformation of the 
wavefunction which can turn it to be periodic, and so the periodic versus antiperiodic boundary condition can be 
used for a topological classification (Z^). Note also that once G(4>) is given between and tt, its value between tt and 
2tt would then be given. 

For our specific example, we can work out easily 



and so 

ci(A\ = ( e o 

where r 3 = ^ q ^ ' so ^K^) ^ s an tip er iodic when — > + 7r. We reproduce the well-known result that our p-wave 

planar state is Z2 non-trivial. 

Our result for G(<fi) can be easily understood, since the two spin species are decoupled. The factors e T1< ^ in its 
diagonal elements just reflect the =Fl Chern numbers for the two spin species. The same calculation for an s-wave 
superconductor would show that G(0) is periodic (in fact a constant), hence topologically trivial. This example, 
though elementary, would be useful to understand the more complicated cases below. 

B. Superconductor with Rashba interaction 

We now consider a superconducting system without up-down reflection symmetry This system has now been 
studied much theoretically (e.g i 26 ' 27 i 28 ' 29 i 40 ' 41 ). First let us consider the normal state. In addition to the kinetic 
energy Hk quadratic in k, we now add the Rashba (single-body) term 

H R = ^2 - a R{ k )a\ a {z x k) ■ aa.pa^p (18) 

Here a are the Pauli spin matrices, and otn(k) > is a measure of the strength of this Rashba "interaction". This 
term can be regarded as a k dependent magnetic field, pointing along z x k, thus the system is still rotationally 

symmetric about its normal. We shall take an(k) — > as k — > 0, so that Hr is well-defined there. In the normal 

,2 

state, the eigen- energies can easily be seen to be £± (k) — =p aji(k) — fj,, according to whether the spin is parallel or 

anti-parallel to the effective field at k, and we shall denote these two bands as ±. These two bands thus have different 
Fermi radii kp±. The spin dependent part of the wavefunctions can be chosen as 



> 



(19) 

1 / ie'^ 
V2\ 1 



We note that 6|fc+ >= ie" l ^| - k+ > and Gljfc- >= -je^sl - k- > and thus 9 at 9" 1 = ie^a 1 - and 

11 11 k+ —k+ 

9ot 9 _1 = —ie^^ka} - . We note that |fc± > do not form a good pseudospin basis since they do not satisfy eq ([5]) 

and ©. We shall refer ± as helicity indices. 

For the superconducting state, in the limit where the pairing is weak compared with the splitting between the two 
Fermi surfaces, it is sufficient to consider only pairing terms on the same band 4 ^. Thus we can write 



-^pair 



\ J2 (iA + (k)e~^al + al u - iA_(fc)e^^4_a t _-_ + c.c.) (20) 
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Huau- is time reversal invariant if both A±(fc) are real, since we have, e.g., Oie ^ial a* _ 1 = ie a} - . 

pau ^ K ' ' ' 6 ' fe+ -fc+ k+ -k+ 

The origin of the e^ % v% terms here should not be confused with those in eq ([2]): they are the result of the choice eq 
(|19p and do not necessarily imply p-wave pairin g 26 ' 41 . Using eq (fT9")) we can rewrite H pa i r in terms of ordinary spins 
and find 



(21) 



with 



A(k) 



^ A + (fc)-A-(fc) c _j, 

A + 2 (fc)+A-(fc) 
2 



A + (fc)+A_(fc) 



;A + (fc)-£. 



(22) 



Thus it is a linear combination of s and p-wave planar pairing, with A s (fc) = ^+iM+^ ( fc ) anc j A p (fc) = A +W 2 A — . 
If A + (fc) = A_(fc) for all k, we have pure s-wave pairing, and if A+(fc) = — A_(fc), i? pa ir reduces to i^i""- Generally 
A±(fc) = A s (fc)± Ap(fc). For ease of discussions below, we shall always assume that A+(fc) is positive, where as A_(fc) 
may change sign as a function of fc. 

As k — ► oo, all pairing terms are negligible compared with the kinetic energy k 2 /2m, and so the system is effectively 
normal. To avoid the complications as mentioned below eq (|19p we shall require further that an(k) — > there so that 
we can easily write down a good pseudospin pair 



DX 



(1,0,0,0) 



(23) 



vff 1 = (0,1,0,0) 



0, since A p vanishes, 



which is in fact the same as |\I> A > defined in the last subsection in this k — > oo limit. At k 
the pairing either vanishes completely, or is always s-wave like only. Since ctn(k) there vanishes as well, we can have 
a pair of good pseudospin wavefunctions, which we write as 



*f 2 = N D2 



D2 



N 



D2 





-A,(0)/(25(0)-£(0)) 
1 




/A s (0)/(i?(0)-£(0)) 





(24) 



V 



w ith N D2 = (( E(0) - i{Q))/2E{Q)) 1 /' 2 . Here £(0) = £±(0) = -ft is the kinetic energy as k -> 0, and E(0) = 

For general k, the Hamiltonian can be easily diagonalized in the ± helicity basis. The energies are ±E±(k) with 
E±(k) = \J^±(k) + A 2 _(fe). The wavefunctions in ordinary spin and particle hole space can then be found by the 
transformation eq (fT9|) . We can choose 



\ 



-ie**iA + (k)/(E+(k)+£ + (k)) 

V A + (fc)/(s + (fc)+e+(*)) / 



(25) 



and 



N 



ci 



V2 



( 



ie 



\ 



-A_(*)/(E- (*)+£-(*)) 



(26) 
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where N^ 1 (k) — ((E±(k) + £,±{k))/2E±{k)) 1 ^ 2 are normalization constants. These functions are chosen so that they 
reduce to the normal state wavefunctions eq (| 19[) when k — -> oo. They are well-defined for finite k (0 < k < oo) if A_ (k) 
does not change sign anywhere for k < kp-. (If A_(fc) changes sign at fc* < &f-, since A-(k)/(E-(k) + £_(&)) = 
(-EL(fe) — £_(fc))/A_(&) and < 0, {E_(k) - £_(fc))/A_(fc) — > ±oo according to whether fe approaches fc* from 

the sgnA_(fc) > or < side.). However, there is no harm if A_(fc) changes sign for k > kp-. 
There is an alternative choice of wavefunctions 



qC2 



(k) 



JVf 2 
V2 



( A + (k)/(E + (k)-Z+(k)) \ 



V 



(27) 



/ 



and 



^ 2 (k) 



N 



C2 



V2 



( ie-**sA_(jfc)/(£_(*) - £_(*)) \ 
A_(fc)/(£_(fc)+^_(fc)) 
-1 



V 



(28) 



where N^ 2 (k) = ((E±(k) - £±{k))/2E±(k)) 1 . These wavefunctions are well-defined for < k < oo only if there 
are no sign changes of A_(fc) for k > kp-- (Actually > and |\&2 2 > are always identical, since we assumed 

A+(fe) > 0. We just write explicitly j'1'9 2 > again for symmetric purposes.) 

To study the topological classification as in subsection III A[ we need to construct good pseudospin bases satisfying 
eq (JS|) and © for different parts of k space, and then try to match them at the boundaries. However, Ivf^ 1 '* 72 > 
do not form good pseudospin bases. The problem is that, as in the normal state, the time reversed partner of the 
+ band is the + band itself. Therefore to have a good pseudospin basis one must have \^i(k) > proportional to 
|4 , + (fc) > over part of the k space whereas it is proportional to |^_(fc) > in the other, introducing discontinuities in the 
wavefunctions. Rather than trying to deal with this rather messy situation, we note that the topological classification 
must be independent of the specific construction of the pseudospin basis. If one has a smoothly defined wavefunction, 
such as l^j- 1 ' 2 > above under suitable circumstances, then we should be able to decide the Zi invariant from these 
wavefunctions alone. We propose that this can done by studying the SU (2) part of the matrix 



T(4>) =< * D1 (<^)|* C (0) > fe ^ 00 < * C (</))|* D2 (0) > 



fc^O 



(29) 



with the matrix elements of each matrix defined as in eq ©■ Here 1 — > +, 2 — > — for \*f?± >, and C — > CI or C2 
should be chosen so that it is smooth for < k < oo. Intuitively, we are trying to consider the transformation between 
the good pseudospin bases |\& D1 > at k — ► oo and |\I> D2 > at k — * through intermediate states |<P± > which are 
smoothly defined for < k < oo. Eq (j2U|) is a special case of a more general formula (|Aip which we shall explain in 
the Appendix. 

so l*^ 1 > 



(30) 



We now evaluate eq ([29]) . Let us first consider the case where A_(fc) never changes sign for k < kp^ 
can be used in eq (f2"9"| . We obtain 



( 1 



-i<p 



independent of the sign of A±(fc). For k — > 0, we get 

< ^2,^1 0= ±(~ie i * -sgnA_(^_) \ 



(31) 



Here we have used the fact that sgnA_(fc — > 0) = sgnA_(fcp_) since A_(fc) never changes sign for k < kp-. Thus, 
for sgnA_(fci?_) > 0, we get 



T(0) 



2 V ie 1 ? 1 

1 

-1 



1 

-1 -ie l 



(32) 
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independent of <j>. This state is Z2 trivial. For sgnA_(fcp _) < 0, we get instead 

ru\ 1 ( 1 ze " # ^ f K ~ l4 1 
A(0j _ 2\ie* 1 J { 1 ie* 

ie " 10 ^ f33) 

the same as eq (JTHJ) , the result for the planar phase in the last subsection. This state is therefore Z2 non-trivial. We 
note that the above calculation applies even when A_(fc) changes sign for some k > hp-- 

For the case where A_(fc) changes sign for k < kp- but not for k > kp-, we must employ instead \*&± 2 >. We 
find 



1 >k ^°°- y/2\ie* sgnA_(fc,_) 
where we have used sgnA_(fcf_) = sgnA_(fc — * 00). On the other hand, 

,D2„ TI C2. 1 f-ie*+ -1 



(34) 



independent of the sign of A_(fc^_). Hence, if sgnA_(fc F _) > 0, we get 

™ = 2 U # 1 A -1 - 
1 
-1 

as in eq (|32|) . independent of 0. This state is Z2 trivial. If sgnA_(fcp_) < instead, we get 

1 ( 1 -ie-t* \ ( ie~ l * 1 
±W _ 2 V*e** -1 J\-l -ie* 

ie-W 
ie* 



(36) 



(37) 



as in eq (j3"T|) . This state is again Z2 non-trivial. 

We thus conclude that the system is Z 2 trivial or non-trivial according to sgnA_(fcp_) > (<)0, correspond- 
ingly whether A s (k) > (<)A p (fc) at kp-, a criterion which we had obtained before by considering the vortex 
bound states^. This result can also be understood from continuity. When changing A-(kp-), from say, posi- 
tive to negative, one must go through the point where A_(fcp_) = 0, which indicates a gapless phase since then 



E-(kp-) = \l £L(kp-) + A 2 _(kp-) would vanish on the Fermi surface £_(fci?_) = 2 -. Hence if A^(kp + ) and 

A_(fc F _) are of the same sign, then one can smoothly change A±(k) to obtain the s-wave state A+(fc) = A_(fc) 
without going through any gapless phase. Similarly, if A-(kp+) and A-(kp-) are of the opposite sign, then the state 
is connected to planar state A + (k) = — A_(fc). 

Fu and Kane^ also proposed that studying the matrix elements (more precisely, the Pfaffian) of between 
wavefunctions at the same k point can distinguish whether the system is Z2 trivial or non-trivial. For the system 
with only two (un-)occupied states \^i y 2(k) >, this just amounts to studying the function P(k) =< ^2{k)\Q\^i(k) >■ 
The non-trivial topology of the system is determined by the presence of odd numbers of pairs of point zeros of P{k). 
Since this overlap is evaluated at a given k (not relating ±fc) and a transformation of basis would not affect the 
presence or absence of zeros of P(k), ^1,2$) > here need not be well-defined pseudospin states in the sense of eq (O 
and ([6|). From our previous calculations, we can also see in our example how P(k) can distinguish Z2 trivial versus 
non-trivial cases, provided one can use a smooth wavefunction in k for l^i^fc) > (Here |4'^ 1 ' C ' 2 (fc) >), except here 
that zeros of P{k) must now appear as circles. Suppose that A_(fc) never changes sign for k < kF~, then we can 
employ |*± X (A) >. We find 

P(k) = < * cl (fc)|6|*5 1 (fc) > 

_ N Cl N Cl r-i , A +( fc ) A -( fc ) 1 /oox 

" N+ ^ [1+ E + (k) +M k)E^k)+Uk) ] (38) 
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As k -> oo, P(k) = 1, but as fc -> 0, P(fe) -> sgnA + (0)sgnA_(0) = sgnA_(0) = sgnA_ (fc f _) 42 . If A_(fc F _) > 0, 
then .P(O) = 1. However, if A_(/cp_) < 0, then P(0) = — 1, implying P(k) must vanish at some intermediate k 
since > are smooth. In the special case /cf+ = kp- = kp, this zero occurs exactly at /cf, since there 

£+(fc) = £,-(k) = and so P(kp) = [1 + sgnA + (fcF)sgnA_ (kp)}/2. If fc^+ and kp- are not equal, it can shown that 
P{k) vanishes somewhere between these two Fermi surfaces. 

Our calculations above do not apply if A_(fc) changes sign for both k > and < kp-- If this happens, a further 
generalization of eq ((29)) is necessary. We can now write 

T(0) =< * D1 (0)|*C1 (0) >oo< *C1 ( ^|^C2 (0) >feF _< *C2 (0) |^D2 W >q (3g) 

This equation is a special case of eq (|A13|) in the Appendix. Using eq (|5D|) , (f3"5)) and 

<* C1 WI* C2 W>^=(; sg nA^_)) ( 4 °) 

we can easily check that the SU(2) part of T is periodic or antiperiodic in <f) — > <fi + n depending on whether 
sgnA^(fc f _) > (<)0. 

C. General spin-orbit interaction 

We generalize the above to general spin-orbit interaction 

H so = ~a\ a h ■ ^ a ,fiau,p ( 41 ) 

We parameterize the direction h(k) by the spherical angles (/3r, ag), and thus we have hr — 
\hj:\(siii[3j:cosaj:, sin/J^sinag, cos/3g) in Cartesian coordinates. Time reversal invariance requires 

\h(<f> + w,k)\ = \h(<l>,k)\ 

0(</> + vr, k) = 7r — f3((f>, k) (42) 
a(cj) + n, k) = a((f), k)+Tt (mod 2tt) 

but we shall otherwise allow general fc dependence for these quantities (cylindrical symmetry would not be enforced). 
\h^\ will be required to vanish at k = 0, and we shall assume that it is finite for < k < oo so that the bands are 

non-degenerate. To have well-defined a^, /?£ cannot be equal to or 7r. We shall assume that hj: does not cover the 
entire sphere, so with suitable rotations of the spin quantization axes, this condition can be fulfilled. Otherwise some 
extra care has to be exercised when using the formulas below, but we shall not go into these complications here. If 
\h^\ is independent of directions, and further (3^ = ir/2, — 4>j: + ir/2, then (|4"Tj) reduces to the Rashba interaction 

d. 

For the normal state, the eigenfunctions, with energies £±(k) = =p \h~\ — fj,, can be chosen to be 



|*+ > 



cos(/%/2) 
sin(/%/2)e M £ 



\k-> = [ -<**Wi/2>- ia ' 



(43) 

cos(/%/2) ) 

generalizing eq (fl9|) . We now have the time-reversal relations 0|fc+ >= — e~ la k | — k+ > and Q\k— >= ~e la k | — k— >. 
The angular dependent Fermi surfaces will be denoted as kp±(cf>). 

In weak pairing limit so that the pairing is only between degenerate bands, we have 



-ffpair , } 



\ £ (-A + (fc)e--%+ + a t _ s+ - A_(fc)e^4y_._ + c.c) (44) 



k 
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generalizing eq (|20|) . Here A±(fc) are real. As before, the e Tia x phase factors are due to the wavefunctions |fc± >. 
Hpair obeys time-reversal invariance. It is a mixture of singlet and triplet pairs, but the later is not necessarily a 
simple planar state. We shall take A + (k) > always, but A_(k) can take either sign. If the system is to remain 
gapped, then A_(&f_(0)) must be either be all positive or all negative for all angles 0, a situation which we shall 
assume. Note that as k — > 0, all components of the p-wave pairing must go to zero, with only the s-wave component 
A s surviving. There, A+(0) = A_(0) = A s (0). 



In the limit k — > oo (0), we again have the good pseudospin basis |\I> > (I* >) as in eq (j2"5)) (eq For 
< k < oo we can solve the Hamiltonian again first in the ± helicity basis, then perform the transformation using 
(l43l). These wavefunctions can be written as 



*% 1 (k)=N 



ci 



/ cos(/%/2) \ 

sin(/%/2)e iQ s 
-sin(/3 fe -/2)e-cA + (fc)/( J B + (fc) + £+(£)) 
\ co S (/%/2)A + (fc)/(i? + (fc) +£+(£)) J 



(45) 



and 



-sin(/3g/2)e- M £ 
cos(/%/2) 

-cos(/%/2)A_(£)/(£_(£) + £_(£)) 
V -si n (^/2)e--*A_(fc)/(B_(fe) + e-(fe)) y 



(46) 



These wavefunctions are smooth in k if A_(fc) do not change sign for k inside the kF- Fermi surface (fi > 0). 
The alternative choice of wavefunctions 



^% 2 {k) = N 



C2 



( cos(^/2)A + (fc)/( J B + (fc)-e+(fc)) \ 
sin(/%/2) e ^A + (£)/(£ + (fc) - £+(£)) 
-sin(/%/2)e iQ S 
\ cos(^/2) / 



(47) 



and 



C2 



( -sin(^/2) e - M £A_(fc)/( J B_(fc) - Uk)) \ 
cosC9jj/2)A_(fc)/(B_(fc)-e_(fc)) 
-co S (/%/2) 

\ -sin(/3 £ /2)e- M £ / 



(48) 



are smooth if A_(fc) does not change signs for k outside the kp- Fermi surface. 
To study the topology, we consider the matrix 



T(0) ^< * D1 (0)|^(0) > k 



c, 



-iAa(4>)T 3 /2 



(49) 



again as the transformation between | , 3/ D1 > and |SI> D2 > via smooth wavefunctions (C — * CI or C2) defined in the 
intermediate k region. The factor e -^Wr 3 /2 is to take 

care of the variation of the wavefunction \*& c (cf)) > from 
k — > to k — » oo, here Aa{(j>) = a(4>, oo) — a(0, 0). For the justification of this formula, see the Appendix. 
Going through algebra similar to the last subsection, we find, if A_ (kp- (<j>)) > 0, 



TWO 



s in( A S Lzfti) e - i ( Q =o+ a o)/2 cos( /j °°"^° )c"' Aa/2 



— cos( 



/3qo-/3q -| c iAq/2 



sin(i 



■)gi(aoo+Qo)/2 



(50) 



Here /3 C 

interaction, eq ([50)) reduces to 



oo) etc, and we have suppressed the <j> labels in eq (|50|) to shorten notations. For the Rashba spin-orbit 
^1 J We ^ ave me ^ before. For here, we note that detT = 1, so that the SU(2) 

part of T can be simply taken as T itself. Generally, with eq ([42]) , we see that T(0) is periodic when <fi — > <f> + tt, 
as Aa(0 + 7r) = Aa(4>) (see Appendix), and both sin( ;3 °°~^° ) and e ± i ( a =o+ Q o)/2 cnan g e s ig n . This state is therefore 
topologically trivial. 

If A_ (fci?- (</>)) < 0, we find instead 
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/- sin (£oc±&) e -i(a»+ao)/2 cos (£°o±£o) e -iAa/2 \ 
A ^-( v cos (^+& )e iAa/2 sin( ^+4 )e ^ +Q0 )/2 J 

If /3g = 7r/2 and ~ (ft + tt/2, then this formula reduces to ie~ l ^ Ts we have met before. Generally, we note now 
detT(0) = —1 independent of 0, so we can take for example G((ft) = — iT(0). Using eq PS|) , we find that T(0), hence 
G(0), is antiperiodic when (ft (ft + tt, as cos(^ 2£ j^ a ) and e ± ^ Qoo+Q °)/ 2 both change sign, showing that again the 
state is topologically non-trivial. 

Let us also examine P(k) =< <& c (k)\G\V°(k) > with C -> CI or C2. Again |P(fc)| = 1 at jfe = or oo. If 
A_(fci?_) < for all 0, we can show that again P(k) must vanish somewhere in k space. Explicit calculations show 
that P(k) is real, so it still must vanish on a line in k space. For kF+(4>) = kF—{<f>) = kp{(ft), this happens exactly on 
the Fermi surface (which in general is not a circle). 



III. CONCLUSION 



To conclude, we have considered the topology of two-dimensional time-reversal symmetric gapped superconductors 
by directly analyzing the quasiparticle wavefunctions, specifically the transformation function between wavefunctions 
defined in different k space regions. This transformation function allows us to classify these superconductors into 
two types. For weak-coupling superconductors, we verified that the relative sign of the order parameter on the Fermi 
surfaces is the crucial parameter which determines the topology, providing a deeper understanding of earlier results 
on vortex bound states^! and edge states 2 ^ in these systems. 
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APPENDIX A 

Here we explain the formulas eq (|29[). (|49|) used in text. Suppose that we have a well-defined pseudospin basis 
|* (fc) > for k > fci, and another |\I rD2 (fc) > for k 2 > k > 0, here fci > &2- Suppose further that we have a solution 
(but not necessarily satisfying eq (O and ©) |\& c (fc) > to the Hamiltonian for k\ > k > fe. Then we consider the 
matrix 

T(^k u k 2 ) =< * D1 (0)|VI>C (0) >fel e -iAa(*)Va < *C (0) |^D2 (0) >k2 (M) 

Here Aa = a(<f>, fci) — a(<fi, fe), where a(0, k) was defined in section III CI The matrix T is in the same spirit as eq 
© except that we now examine the transformation between the two pseudospin bases I'S' 01 > and |\I> D2 > via an 
intermediate |* c >. The factor e ~ jAQ ( < ^ r3 / 2 is to compensate for possible variations in between k = ki and k — k 2 . 
In this formula, we have limited ourselves to the case where 0% never equals to or tt. This guarantees that the factor 
e -iAa(0)T 3 /2 j g well-defined despite that a is defined only up to mod 2tt: For a given k, if say a((ft, k) increases by tt 
when (ft advances by tt, then the same thing must happen for a neighboring k! since a((ft, k) is smooth (instead of, e.g., 
a{(ft, k') decreasing by tt). Continuing this argument shows that a(<p + tt, k\) — a((ft, k%) = a{(ft + tt, k 2 ) — a((ft, k 2 ), so 
that Aa(4> + tt) = Aa(<p) and Aa(<ft + 2tt) = Aa(<ft) without any 2tt ambiguities. The precise form of this factor is 
dependent on the behavior of the hclicity basis |\f c > under time reversal (see eq (|A5[) below). A more complicated 

form is necessary for example if we choose to multiply eq (|45l) - (|48|) by k dependent phase factors, which we choose 
not to do here. 

For simplicity, we have also written eq (|A1[) on two circles ki and k 2 in k space. Since the system may not be 
cylindrically symmetric, one can also choose to write this equation on more general paths in k space enclosing the 
orgin, but this would complicate the notation. In the text we have only used this formula for k\ — > oo and k 2 — * 0. 

T((f>, ki, k%) is obviously unitary. Let us first write 

t D1 - c (0) = < * D1 (0)|* C (0) > fel 

= e^ m ' C (0)G D1 ^(0) (A2) 



12 



where detG D1,c = 1 etc. We thus have 

01,C ,,-,,C,D2 



T(^) = ^(x iJ ^(0)+ X ^^(0)) G Dl,C ( ^ e - l Aa(0)r 3 /2 G C,D2 (0) 

= e ix ^G(<f>) (A3) 
Let us now investigate the consequence of time-reversal symmetry. Eq (J5J and |6j) require 

fltff^ + T) >,|* 2 m (0 + ^) >) fc = (levi/f 1 ^) >,|e*?V) >) fc ( J) (A4) 

For our hclicity basis, we have, from eq (|45 p -(|48 p 

(1^(0 + 7r)>, + 7r)>) fc = (|e^(0)>,|e* c (^)>)^ e Q _ e _i Wl *)J (A5) 



We thus obtain 



Similarly, 



' c - M (* + - ( -7*' _ e ,i fe) ) <*»»<*))• ( _°, ; ) (at) 

By judicially inserting 1 = T^qJ^^ ^ arL< ^ using the commutation properties of Pauli matrices, we find 

T(4> + 7T) = e -Hx D1 ' C (0)+X C ' D2 (0)) G Dl,C ((/))e - t A Q (0)r3/2 G C,D2 ((/)) (Ag) 

where we have used Aa(0 + 7r) = Aa(</>). We then obtain 

e *x(*+*)G(^ + tt) = e~ lx ^G(0) (A9) 

thus allowing classifications by the behavior of G(0) under — > + 7r as in section Hi Al 

Eq (|A1[) is directly related to the transformation used in£i2I and formula such as eq © . The analogy to these would 
be to study 

t Dl,D2 ( ^ =< ^Dl (0)|,JrD2 (( ^ >R (A10) 

at some K where both |* D1 (fc) > and |* D2 (fc) > are well-defined. If indeed |* D2 (fc) > in eq (|A1|) is defined for k 
up to fci, then we see that, using that |\& D2 (fc) > forms a complete set for the positive energy levels, 

T(0) =< * D1 (0)|* D2 (0) > kl M(<t>,k 1 ,k 2 ) (All) 

where 

M(0,fc!,fc 2 ) =< M/ D2 (0)|VI> C (0) > kl e - iAa ^ 2 < * C (0)|* D2 (0) > k2 (A12) 

so amounts evaluating (|A10|) at ki and post- multiplying by the matrix M. By the same analysis as from eq (IA1[) to 
(|A9[) . we see that the SU(2) part of M must be either periodic or antiperiodic when <fi — > cf> + tt. However, since the 
wavefunctions entering M are smooth and M(</>, k 2 , k 2 ) is necessarily simply 1, the SU(2) part of M((/>, k\, k 2 ) must 
be periodic when <f> — *■ <f> + tt. Hence the classification is not affected by the post- multiplication by this matrix. 

Similarly, one can show that the classification by the SU(2) part of (|A1|) is independent of the choice of the 
pseudospin basis |\J/ D1 > or |\I> D2 >, so long as they are smoothly defined in their respected regions in k space. A 
different choice would just amount to pre- or post-multiplying eq (|A1[) by a matrix with the SU(2) part of which being 
periodic in <f> — > + tt. 

We have assumed that /3g is never or tt in eq (|Aip . When this does not apply, the factor e - lAa W T 3/2 there has 
to be rewritten. The necessary form can be obtained by considering more carefully how eq (| A5|) changes when k is 
varied, but we shall not go into these complications here. 
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When there are sign changes of A_(&) for both outside and inside one must use \^ ci ((f>) > and \*& C2 ((f>) > 
for these two regions separately. Similar considerations above show that a transformation matrix which can be used 
is 

T(&*i,fc) = < * D1 (0)|* C1 (0) > kl e -<A«x(*Wa < *ci W |^c2 (0) >fc ^ w) 

xe - 2 Aa 2 (0)r 3 /2 < ^C2^|^D2^ > ^ (A13) 

where Ai(0) = a(</>, fci) — a(</>, fc^_ (</>)) and ^{(j)) = a{4>, hp— (<fi)) — ct((fi, ^2). If a is independent of /c, this formula 
reduces to eq (|39f at the end of section III Bl 
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(A14) 



Note </>_g = 0^ + 71 mod 2ty. 



When discussing the topology, we regard k — > 00 as a single point irrespective of its direction, as done in-^ etc. We shall 
therfore require well-defined pseudospin wavefunctions to have <f> independent values as k — > 00. This condition would not 
be imposed upon on wavefunctions such as > in subsections III Bl and III CI that are used as intermediate wavefunctions. 
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